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Abstract 

A method-of-moments scheme is invoked to compute the asymptotic, 
long-time mean (or composite) velocity and dispersivity (effective dif- 
fusivity) of a two-state particle undergoing one-dimensional convective- 
diffusive motion accompanied by a reversible linear transition ("chemical 
reaction" or "change in phase") between these states. The instantaneous 
state-specific particle velocity is assumed to depend only upon the instan- 
taneous state of the particle, and the transition between states is assumed 
to be governed by spatially-independent, first-order kinetics. Remarkably, 
even in the absence of molecular diffusion, the average transport of the 
"composite" particle exhibits gaussian diffusive behavior in the long-time 
limit, owing to the effectively stochastic nature of the overall transport 
phenomena induced by the interstate transition. The asymptotic results 
obtained are compared with numerical computations. 

Keywords: Homogenization, Brownian motion, macrotransport theory, 
generalized Taylor dispersion 
PACS: 05.40.Jc, 87.10. +e 

1 Introduction 

Two-state models are often invoked to construct simple, easily analyzed trans- 
port models of otherwise overwhelmingly complex processes. Examples of such 
processes include chromatography ^E] , isomeric conversion of proteins in solu- 
tion and in vivo calcium kinetics and ratchet-driven molecular motor pro- 
cesses pi]. In such models, it is typically assumed that the species of interest, 
either a single entity (dispersed in a passive solvent) or a continuum solute con- 
centration field composed of a dilute collection of such non-interacting species 
(again dispersed in a passive solvent), undergo(es) a microscale, state-specific, 
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convective-diffusive transport process wherein the transport coefficients depend 
at time t upon the "state" (a, b) of the species at that time. The rate of inter- 
change between states is assumed to be governed by first-order kinetics. Since 
the total species number density does not change in time, this kinetic process is 
equivalent in its consequences to a reversible first-order reaction or phase tran- 
sition, a <-» b, between states a and b. In almost all cases of interest, one is not 
generally concerned with the detailed intrastate transport of the species in state 
a or b individually, but rather only with knowledge of the long-time, asymptotic 
aggregate transport of the composite a-b species; that is, interest focuses only 
on the combined transport occurring in both states a and b for times sufficiently 
long such that the interstate reaction has effectively achieved global equilibrium 
(although the reaction may not be in local equilibrium at each point in space). 
Remarkably, this averaged aggregate transport is diffusive, even in the absence 
of molecular diffusion. 

Many two-state models only feature the explicit incorporation of convec- 
tion and a "reaction," without explicitly accounting for intrastate molecular 
diffusion. Consequently, it is tempting to say that these models constitute ex- 
amples of "convective dispersion without molecular diffusion." However, the 
interstate transfer reaction often corresponds to a simple model for the solute's 
molecular diffusion, allowing it to sample two different velocity fields prevailing 
at different points in space, say, in two different phases. In an exactly-posed 
model 0, molecular diffusion transports the solute to the interphase bound- 
ary, whereupon an interstate mass transfer "reaction" occurs. Without the 
explicit incorporation of diffusion, a simplified (or rate limited) one-dimensional 
model |1I2| allows the solute to "jump" between phases in a kinetic manner. The 
dispersivity computed from such a "nondiffusive" model 0, which possesses a 
mathematical structure similar to the generic problem to be considered here, 
does indeed furnish a convective dispersivity which is explicitly independent of 
the molecular diffusivity. From a physical viewpoint, though, the convective 
dispersion may still be attributed to molecular diffusion, but at a finer-scale 
level of description than is represented by the coarse-scale transport equations 
adopted in our model. As such, it is not entirely correct to say that existing 
two-state chromatography models satisfy our claim of "convective dispersion 
without molecular diffusion." 

Rather, our claim more closely corresponds to the case where the disper- 
sion arising from the interstate reaction does not occur as a consequence of 
molecular diffusion involving solute molecules which randomly shuttle back and 
forth between two different velocity fields (phases) located at different points 
in space. To see how such a scenario arises, consider the examples of molecu- 
lar motors and denaturing of double-stranded DNA. In the case of molecular 
motors 0], the species can exist in either a charged or uncharged state. When 
charged, the particle undergoes convective transport due to its interaction with 
a charged potential field, the latter referred to as the "track." During the course 
of its convective motion, the particle may also interact biophysically with the 
surrounding medium (say, by incorporating ATP). The latter interactions neu- 
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tralize the particle and arrest its convective motion. 1 The overall biophysical 
process is cyclical, so the particle may again become uncharged, enabling the 
transition between charged and uncharged states to be modeled as a reversible 
chemical reaction. Though the particle itself also undergoes Brownian motion, 
the transitions between the charged and uncharged states depend only upon 
the chemical interactions of the particle with the surrounding medium, not the 
spatial location of the particle. 

Another example of a non-diffusive change in state involves the partial de- 
naturing of double-stranded DNA. In such circumstances, the two associated 
strands become locally separated over a short range of base pairs, owing to an 
elevated temperature (or an equivalent concentration of denaturant) exceeding 
the local melting temperature. The native-denatured transition is primarily 
attributed to a chemical interaction between the base pairs (and the denatu- 
rant, if present), rather than to molecular diffusion. The native and denatured 
states possess distinctly different electrophoretic mobilities, and it is these dif- 
ferences which have been exploited in the context of numerous electrophoretic 
schemes |7||S] for the detection of disease-causing mutations. Inasmuch as the 
mutation detection relies upon interpreting the chromatographic bands, it would 
be useful to have a greater understanding of the role of fluctuations between the 
states (which could be related to, say, fluctuations in the effective temperature 
of the medium) upon the overall chromatographic process. In both the motor 
and denatured DNA cases, the change in the solute state, and thus the so- 
lute velocity, arises from a chemical interaction with the medium, not from a 
conventional molecular diffusion process transporting the solute to a new point 
in space. Consequently, the convective dispersion process arising here would 
persist in the absence of molecular diffusion. 

As a prelude to a much more general analysis, the main point of this paper 
is illustrated by considering the following elementary non-diffusive system: 

BP BP 1 

~BT + U ~d£ + k ^- Pfc ) = 2 6 {X) 6 {t) ' (L1) 

~bt ~ u ~d + k ^- p ^ = 2 s {x) 5 {t) ' (L2) 

where P a (x, t) is the conditional probability density of state a, U is the species 
velocity, k is the interstate reaction rate, and 5 is the Dirac delta function. Let 
Pa (<7; s ) denote the Fourier-Laplace transform of P a (x,t), with q the Fourier 
variable and s the Laplace variable. Transforming eqs. I|l.ll) - (|1.2|) m t° Fourier- 
Laplace space furnishes the following coupled set of algebraic equations: 

sP a + (iqU) P a + k(P a - P b ) = ^, sP b -(iqU)P b + k(P b -P a ) = ^. (1.3) 

In more realistic models, the tracks possess a spatially-periodic, ratchet-like charge dis- 
tribution. While modeling periodic systems necessitates a more elaborate moment-matching 
scheme than the one employed here, the existing macrotransport scheme [B] for single-state 
systems is readily extensible to two-states. The resulting calculations are more involved, but 
the final result reveals, as would be expected, that a similar contribution to the convective 
dispersion exists in the periodic case. 
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Upon solving for P a and Pb and forming their sum, the Fourier-Laplace trans- 
form of the total probability density adopts the form 

P = P a + P b =(s+^-) . (1.4) 



For long times t 3> 2k, corresponding here to s <C 2k, this becomes 



P«( S + ^) . .1.5) 



2k 

The latter expression is readily inverted to yield the gaussian distribution 



P w (47rL»*t) 1/2 exp 



a; 2 1 



4D*i 

with an effective diffusion (dispersion) coefficient 



(1.6) 



U 2 

D * = W M 

Moreover, the mean position of the composite particle remains at the initial 
position x = 0. The latter property is equivalent in its consequences to a zero 
composite particle velocity, U* = [compare eq. (|3.5H ]. 

In conventional long-time average "macrotransport" analyses |S] of this type, 
the dispersivity D* arises via the stochastic sampling of local velocity inhomo- 
geneities by solute diffusion in the continuous local space. For example, in the 
classic case of Taylor dispersion ;9, , quantifying the global solute transport oc- 
curring in a Poiseuille flow within a long circular cylindrical tube, radial diffusion 
across the streamlines enables a solute molecule being transported axially by the 
Poiseuille flow to sample the fluid's transversely parabolic axial velocity profile 
innumerable times as it moves downstream, eventually attaining a stationary 
mean velocity U* and dispersivity D*. In contrast, the source of the disper- 
sion in our two-state problem resides in the interstate reaction, which allows the 
composite a-b species to sample the discrete "state space." Nevertheless, disper- 
sion arising from the sampling of different state-specific velocities is completely 
analogous to the dispersion which arises from so-called "local-space" velocity in- 
homogeneities, e.g. Poiseuille flow, in classical macrotransport theory 6 . The 
stochastic foundation of chemical reactions is well established P0]> inasmuch as 
the reaction rate k may be interpreted as reflecting the probability of a reaction 
occurring in the interval between times t and t+8t. Moreover, whereas chemical 
reactions are considered to constitute purely deterministic processes in the full 
position-momentum space ( "local space" ) , the process is nevertheless stochastic 
in the reduced space ("global space") employed in conventional kinetic mod- 
els ^T]. Likewise, the comparable reduction from transport in the distinct 
states a and b to transport in the composite a-b state results in a stochastic 
(diffusive) composite transport process. 
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In the present contribution, we investigate more thoroughly this surprising 
property of two-state systems. This is effected by computing, via the the- 
ory of macrotransport processes (Sj, the asymptotic, long-time mean velocity, 
{/*, and dispersivity, D* , of a more general two-state system than that de- 
fined by eqs. i|l.l|) - (|1.2|) . This scheme may be likened to a multiple-time scale 
analysis of the phenomena, where the short-time, two-state behavior ulti- 
mately serves to determine the long-time, state-independent temporal behavior 
of the system as a whole. In a very different, non-reactive context than that 
considered here, moment-matching concepts have been used to consider 

systems where the different "states" correspond to the different physical phases 
(e.g. solid and liquid) through which a solute molecule can be transported at 
different velocities. In contrast with the present analysis, the latter analyses em- 
ployed molecular diffusion to transport the solute to the interface across which 
the interstate mass transfer (change in state) occurs. Subsequently, Iosilevskii 
and Brenner |14j used eigenfunction expansions of the moments to develop a 
general macrotransport scheme for the analysis of reactive mixtures in an in- 
compressible solvent flow, where molecular diffusion again played an important 
role in enabling the solute to sample the local space. Balakotaiah and Chang 2 
used center-manifold techniques to examine problems in chromatography theory 
which are similar to those considered here, albeit with an implicit dependence 
upon molecular diffusion. The results of the latter correspond to a special case 
of the subsequent analysis, which, we believe, employs a more straightforward 
technique than is manifest in center manifold theory. 2 



2 Generalized Problem Statement 

Consider the conditional probability density P a (x,t \ (p a ) (a — a, b) that a non- 
interacting collection of particles exists in state a and are present at position x 
at time t, given their initial impulsive introduction into the unbounded system 
(—00 < x < 00) at x — and time t = 0, with <p a the fraction of the particles 
initially in state a. The two <p a are not independent parameters, since <p a + 
4>b = 1- If an effective or composite a-b "particle" description is to exist, the 
final results for U* and D* quantifying the spatio-temporal transport of this 
composite "particle" must (and will) prove to be independent of the arbitrary 
choice of labels a and b. Moreover, the macrotransport parameters U* and D* 
will further prove to be independent of the particle's initial arbitrary position 
x = 0, as well as of the initial state fractions (p a . 

The spatio-temporal evolution of the conditional probability densities, P a = 
P a (x, 1 1 <p a ) and Pb = Pb {x, 1 1 </>&), are governed by the coupled set of convection- 
diffusion-reaction equations, 

BP BP B 2 P 

~dT + Ua ~^ ~ Da ~d^~ + k (KaPa ~ KpPp) = (2 - x) 



2 It should be noted that the moment- matching scheme employed here is only valid for linear 
(first-order) reactions, whereas the center-manifold technique may be employed to analyze 



nonlinear transport processes as well 
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valid for a = (a, b) (a ^ /3), with S the Dirac delta function. Here, the constants 
U a and D a represent the state-specific velocity and molecular diffusivity of the 
particle. It is assumed that the latter parameters depend solely upon the 
state (a, b) of the particle, being independent of position x, as well as of time 
t. While it is possible in the interest of simplicity to set the state diffusivities 
D a to zero at the outset, we retain these terms in order to assess what will 
prove to be their relatively straightforward impact upon the averaged long-time 
transport process [cf. eqs. (|3.8l) - (|3.9[l ]. 

In the above, the transition between states has been assumed to occur via 
first-order kinetics, with kK a representing the transition rate from state a to 
b, and kKb the transition rate from b to a, with all reaction parameters being 
positive: (k, K a , Kb) > 0. We employ the trio of parameters fc, K a and Kb (of 
which only two are independent, say kK a and K = Kb/K a ), rather than a pair 
of reaction rates, say, k a = kK a and kb — kKb, so as to permit a clear distinction 
to be made between the kinetic rate of interstate transfer, embodied in k, and the 
stationary (equilibrium) partitioning of states, embodied in Kb/K a . Moreover, 
this egalitarian notational choice serves simultaneously to emphasize that the 
final results are independent from the arbitrary choice of labels a and b. 

This coupled set of partial differential equations can be solved, in principle, 
subject to the attenuation of the conditional probability densities and fluxes 3 
at infinity for all times t > 0: 

P a ,J a ^0 as |x| — > oo. (2.2) 

Summing eqs. (|2.1|l over a and b, and subsequently integrating over the interval 
(—00, 00) furnishes the particle conservation relation, 

/oo 
(P„ + P b )dx = l ( t > 0) , (2.3) 
-00 

which reflects the fact that each of the particles originally introduced into the 
system at time t = are conserved for all time, independently of the initial 
distributions, <j) a and <pb- 

3 Macrotransport Analysis 

In Appendix A, we use standard moment-matching techniques |S] to analyze 
the generalized problem for asymptotically long times, namely those satisfying 
the inequality 

t > [fc {K a + Kb)]- 1 . (3.1) 

3 When the D a are non-zero, the P a are expected to decay exponentially fast as \x\ — * oo, 
whence the flux attenuation condition in eq. 12.21 will be automatically satisfied. For the case 
where the D a are identically zero, the governing equation 12.11 is hyperbolic, whereupon the 
probability density is identically zero for all times as \x\ — ► oo. 
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Explicitly, the macrotransport parameters U* and D* are computed by match- 
ing the asymptotic rates of change of the total moments of the microscale prob- 
ability density, 

/oo 
x m (P a + P b )dx (m = 0,1,2,...), (3.2) 
-OO 

against the corresponding moments, 

/oo 
x m P (x, t) dx (m = 0,1,2,...), (3.3) 
-oo 

of the macrotransport equation, 

OP - dP - d 2 P 

« a? -'WW- < 3 ' 4 > 

In the latter, P is the so-called macrotransport probability density of the com- 
posite a-b particle, governing the long-time asymptotic evolution of the true 
probability density, P = P a + P b - Solving equation (|3.4|) for P furnishes the 
gaussian distribution, 



P = UnDH) 1/2 exp 



(x-u*ty 



(3.5) 



With respect to the parameters appearing in eq. (|2.1|) . application of the 
moment-matching scheme furnishes the composite Lagrangian velocity, 



T -r* K b U a + K a U b KU a + U b 

u = K a+Kb - (3 - 6) 



where the parameter 



K ^ (3.7) 



constitutes an equilibrium constant, representing the condition eventually achieved 
in eq. I|2.1(l when K a P a = K b P b , i.e. when the "forward" and "backward" rates 
at which the states change are in balance, at least in the global sense embodied 
in eq. (|A.10(I . Consequently, the composite particle velocity H3.6|) represents a 
weighted average of the respective state-specific velocities, U a , with the weight- 
ing corresponding to the stationary partitioning of states at long times. 4 Albeit 
in a very different context, this result agrees with the comparable two-phase 
result of Aris (sans areal factors) [S] , while also reducing to the results of Bal- 
akotaiah and Chang upon setting one of the two state velocities to zero. 
Importantly, the mean velocity depends solely upon the stationary partitioning 
of the particle states, K, rather than being dependent upon the kinetic con- 
stant, k, governing the rate at which the particle transits between states before 
achieving this stationary partitioning of states. 

4 Substitution of the parameters appearing in our initial example 11.11 - 11.21 into eq. 13.61 
confirms that U* = in that case. 



7 



The moment-matching scheme also furnishes the composite particle disper- 
sivity, 

D*=D M + D C , (3.8) 

where 

- M K b D a + K a D b _ KD a + D b 

~ K a + K b = l + K ' (3 ' 9) 
represents the molecular or "Aris" contribution to the total dispersivity, and 

D c = K * K » 3 ^-^ 2 . — 3 (^-^) 2 (3.10) 
[K a + K b f k (l + K) 3 kK a 

constitutes the corresponding convective or "Taylor" dispersion. Of course, 
eqs. H3.8I) - I|3.IU[) reduce appropriately to l|1.7l) in the special state-specific cir- 
cumstances quantifying the latter. It is readily confirmed that D c is positive 
semi-definite, vanishing only when: (i) the state velocities are identical in mag- 
nitude and direction; or (ii) one of the two parameters, K a , vanishes. The latter 
corresponds to the trivial transport process where, for long-times, the particle 
exists in but a single state. The functional form of eq. (|3.1()(l is equivalent 
(to within a prefactor) to that derived by Aris pj- Balakotaiah and Chang's 
result is again recovered by setting one state velocity to zero. 

At the present juncture, it is informative to compare the convective disper- 
sivity, eq. <|3.10(l , with the classical result [§lll5| for the convective contribution 
to the Taylor-Aris dispersion of a solute entrained in a Poiseuille flow moving 
at mean velocity v through a cylindrical tube of radius R, namely 

S° = l** (3.11) 
48 D ' V ; 

where D is the molecular diffusivity of the solute. The quadratic dependence 
upon the velocity "inhomogeneity," (U a — U b ) 2 , appearing in l|3.10|) is analo- 
gous to the functional dependence, v 2 , appearing in (|3.11() . Moreover, the 
factor appearing in H3.I()(I possesses an interpretation identical to the factor 
R 2 1 D appearing in (|3.11|l . Explicitly, each quantity respectively represents 
the characteristic time required for sampling the local space, either the "state" 
space in the present reactive problem or the physical space embodied in the 
tube's cross-sectional area, nR 2 . 

In contrast with U* , which is independent of the kinetic properties embodied 
in k, the convective dispersivity D c depends upon both the particle's equilib- 
rium and kinetic properties, respectively embodied in the parameters K and k. 
The physical basis for this kinetic dependence may be rationalized by considering 
the transport of two particles, each occupying identical positions x and states a 
at time t. Were both particles to remain in state a, deviations in their relative 
mean positions would arise only from molecular diffusion, a phenomenon which 
is accounted for by the molecular contribution, D M , to the dispersivity, D* . In 
contrast, were the first particle to switch to state (3 while the second remained in 
state a, these two particles would tend to separate ("spread") due to the relative 
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difference in their state-specific particle velocities, \U a — U@\. Consider first the 
case where the reaction rate is rapid relative to the rate of convective spreading, 
i.e. k 3> \U a — Up\/Ax, with Ax a characteristic separation distance. [The 
combination fc(Ax) 2 represents the coupling of interstate (local) transport to 
global transport, which is identical to the physical interpretation underlying 
D c .] For fast reactions, the likelihood is great that the particles will return to 
an identical state (either a or /?), and hence identical velocity, in a brief period 
of time, thereby minimizing their "spreading." In contrast, for slow reactions, 
k <C \U a — Up\ /Ax, there exists a high probability that the particles will travel 
significant distances before returning to identical states (and identical veloci- 
ties), whereupon the "spreading" caused by the interstate transfer is expected 
to be much larger than in the fast reaction case. The appearance in D c of 
the square, (U a — Ub) 2 , of the velocity difference rather than the term \U a — Ub\ 
invoked in the preceding scaling arguments, arises from the fact that the disper- 
sivity represents an asymptotic measure of the mean-squared deviation of the 
relative particle positions, rather than simply their absolute deviation. 

An alternative rationalization of the dispersivity D* is achieved by consid- 
ering its role in quantifying the deviations of the instantaneous solute particle 
position, x, at time t from its composite position, U*t, which the particle would 
occupy at time t if it always moved uniformly at the velocity U*. To aid in this 
analysis, define the convected coordinate variable, 

x* d = x - U*t. (3.12) 

Appropriate conversion of the microscale equation p.ljl into this altered coor- 
dinate format furnishes the following equation governing the respective (a, b) 
transport processes: 

BP - BP B 2 P 

^f + ( U <*~ U*) -^-D a -^ + k(K a P a - KpPp) = 4> a 8 (x*) S (t) , (3.13) 

where, here, d/dt = (d/dt) ,. In this new coordinate system, one of the two 
state-specific velocities takes place in the +x* direction, whereas the other oc- 
curs in the — x* direction (except for the trivial case where both state velocities 
are equal), since U* represents an average of the respective state velocities, U a . 
Consequently, the transport process occurring in the convected x* coordinate 
system may be envisioned as a biased random walk, where the state-dependent 
step-sizes are proportional to the velocity difference, \U a — U* , while the prob- 
ability of taking a step in the a-direction is proportional to Kp. 



4 Comparison with Numerical Results 

Numerical solutions are presented here for several illustrative choices of the 
transport and reaction parameters, comparing both the short- and long-time 
evolutions of the exact microscale probability density, P = P a + Pb, with its 
asymptotic counterpart, P, given by eq. 13.51) . For simplicity, we have elimi- 
nated intrastate molecular diffusion effects by setting D a = Db = in (|3.13|) . 
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Anticipating that, for long times, the total solute probability density will be 
convected at the rate U* , as in l|3.6[) . we proceed with the formulation l|3.13[) 
of the microscale transport problem in the convected coordinate system, x* . In 
rendering eq. I|3.13[) dimensionless, it proves convenient to abandon the egali- 
tarian notation employed thus far, choosing state a as the base state. 5 With 
dimensionless time and length variables chosen as 

t = tkK a , x*=x(^), (4.1) 



Ua 

eq. (|3.13|) adopts the respective dimensionless forms 
BP - BP 

~W + Ua W + {Pa ~ KPb) = AlS (r } 6 ® ' (4 ' 2) 

BPu ~ 8 Pi. 

•ir + Ub w ~ {Pa ~ KPb) = (1 ~ M 6 6 ^ ■ (4 - 3) 

The equilibrium constant K is given by (|3.7|) . whereas the dimensionless state- 
specific velocities possess the respective forms 

^ a= lT¥' U b = -KU a , (4.4) 

with U the velocity ratio, 

U d ^ ^. (4.5) 

In the convected, dimensionless coordinate system, our asymptotic result (|3.5J) 
adopts the form 



P = UwD*tj exp 
with D* , the dimensionless dispersivity, 



4D*t 



(4.6) 



k(i-u) 

D* = — ± (4.7) 

{l + Kf 

The coupled equation set (|4.2I) - (|4.3|) was solved numerically for various values 
of the parameters K, U and <f> a , using upwind finite differences for the advected 
terms and forward Euler time integration. Numerical results obtained for short 
it = 0.5) and long (t = 50) times are depicted in Figs. 1 and 2, respectively, along 
with the corresponding asymptotic macrotransport results. Similar results (not 
depicted) were obtained for other choices of the parameters K, U and 4> a . 

The results shown in Fig. 1 demonstrate that, for short-times, the system 
may exhibit dramatically different behavior from that displayed at longer times, 

5 Equivalent results may, of course, be obtained by choosing state b as the base state and 
interchanging all subscripts. 
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depending upon the particular choices made for the parameters K and U and 
initial state 4> a of the system. In general, the short-time transport of a and 
6 as a whole is, not unexpectedly, bimodal, since that species whose veloc- 
ity exceeds U* moves in the +i* direction, while its counterpart moves in the 
—x* direction. Consequently, our asymptotic macrotransport analysis, which is 
symmetric about the origin x* = 0, invariably proves highly inaccurate when at- 
tempting to capture the short-time behavior, typically overestimating the width 
of the distribution. In contrast, as depicted in Fig. 2, the macrotransport de- 
scription agrees very well with the numerical solutions for long-times for all 
choices of the parameters examined. As the species separate, the interstate 
reaction serves to "mix" their probability densities, ultimately giving rise to an 
asymptotically gaussian distribution. 

It should be noted that the hyperbolic structure of the non-diffusive mi- 
croscale equation l|4.2jl - i[4.H|l guarantees that, except for the special case U a — 
Ub, the exact solution of the microscale equations will be asymmetric with re- 
spect to the origin. Explicitly, at time t*, the maximum, ij^axi an< ^ minimum, 
x* nin , possible spatial positions with non-zero probability densities are given 
respectively by ij^ax = U a t* and x^m = Ubi* (taking U > 1 without any loss 
of generality). Aside from the aforementioned special case, the exact solution 
is expected to be asymmetric, since i^ax 7^ ^min- I n contrast, the macro- 
transport solution is not only symmetric about the origin, but also predicts a 
non-zero probability density for all values of i*. However, this asymmetry- 
induced disparity existing between the exact solution and the macrotransport 
solution occurs in the tails of the probability density, at which positions the lat- 
ter density is already exponentially small, and therefore of the order of the error 
incurred in our asymptotic analysis. In contrast, the differences between the 
exact and the macrotransport solutions are expected to be small in the "central 
region" for long times, as confirmed by the numerical results. 
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A Moment-Matching Scheme 

Define the state-specific m*' l -order "local" moment of the probability density of 
state a = (a, b) as 



The latter moments are equivalent to the local space (bounded) variables in con- 
ventional macrotransport theory In order to guarantee that these higher- 




— oo 



x m P a (x,t\cj) a )dx (m = 0,1,2,...). 



(A.1) 
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order moments are finite, it is necessary to strengthen the attenuation con- 
dition (|2.2|l . such that the probability densities and fluxes decay faster than 
algebraically, namely 



P a — > and 



dP a /dx^0 as 



(A.2) 



The differential equation governing the state-specific local moment (|A. II is de- 
rived as follows: (i) multiply eq. (|2.1() by x m ; (ii) integrate over the range 
— oo < x < oo; (iii) integrate the resulting expression by parts; and (iv) sub- 
sequently apply the attenuation conditions (|A.2|) . This scheme eventually fur- 
nishes the relation (for a ^ 0) 



dPj m) 
dt 



--mU a Pi 



m-l) _|_ 

4> a s(t)s m0 , 



(A.3) 



where <5 m o is the Kronecker delta function (i.e. <5oo = 1 and S m Q = for m / 0). 
In principle, this equation can be solved recursively for Pa (t \<p a ), although 
with use of macrotransport theory it proves unnecessary to do so in order to 
eventually calculate U* and D* . 

In a similar manner, the m'^-order "total moment" of the composite particle 
represents the sum of the corresponding local moments [cf. eq. Ij3.2[l ]. 



M m (t | 4> a ) = Pi m) + (m = 0, 1, 2, . . .) 



(A.4) 



Summing eq. (|A.3|I over the two states, a — (a,b), furnishes the following dif- 
ferential equation governing each M m : 



dM n 
dt 



rr p(m-i 
S(t)5 m0 . 



U b P, 



(ro-1) 



+ m(m - 1) \D a P^ + D b Pl m - 2) 



(A.5) 



The macrotransport parameters may then be computed from the following re- 
lationships 



U* = lim 



dMi 
dt 



D* 



lim — (M 2 

{-KXJ dt v 



Ml) 



(A.6) 
(A.7) 



We proceed here to solve eqs. (|A.3|I and (|A.5fl . at least in the long-time limit 
l|3.1fl . This enables us to establish the asymptotic forms of the moments, 
and M m , for m — 0,1,2, thereby permitting U* and D* to be calculated via 
eqs. <|A.6() - (IA.7ll . Upon setting m — in eq. (|A.5|I and integrating with respect 
to t, the zeroth-order total moment is found to be 



M 



(t < 0) . 
(t > 0) . 



(A.. 
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reflecting the conservation condition (|2.3[) . 

The zeroth-order, state-specific local moment equation is obtained by setting 
m = in eq. (|A.3|I . yielding 



+ k (K a P a a - K P P%) = <j> a 5{t) {a =a,b). (A.9) 

The asymptotic solution of this pair of equations, valid for times satisfying the 
inequality (|3.1|> . is given by 

P o«p0,oc +exp) p o,oo = K f (A.10) 

&a + rib 

where the symbol "exp" denotes position- and initial condition-dependent func- 
tions that are attenuated exponentially rapid in time. By forming the sum 

is seen that the state-specific solutions (jA.lOfl 
asymptotically satisfy the normalization condition l|2.3[l . Moreover, for long 
times, the local moments P°'°° are unconditional (rather than conditional) 
probability densities, being independent of both the initial position and the 
instantaneous state of the particle at t = 0, as was to be expected. 

The composite velocity U* is computed by setting m = 1 in cq. I|A.5|) . 
substituting the asymptotic solution (|A.10|) for a — (a, b) into the resulting 
differential equation, and applying eq. I|A.6|) . so as to obtain the result cited in 
eq. (E3J. 

Subject to a posteriori verification, assume the following trial solutions for 
the first-order, state-specific local moments, 

pi M p0,oo (p* t + B ^ + exp; (A n) 

where the B a are time-independent state-specific constants to be determined. 
Form the first moment of eq. (|A.3|I with m = 1 , and use eq. jUD} for a = (a, b) 
together with the trial solution (|A. 1 lfl to obtain 

U* -U a = kK a {Bp - B a ) . (A.12) 

As in prior macrotransport analyses [S] , the latter equation defines each B a only 
to within a common, state-independent arbitrary constant, whose value proves 
irrelevant to the value of D* . Substituting eq. I|3.6|l into (|A.12|) furnishes the 
following relation between the B a : 

The latter time-independent solution, jointly with the fact that summing eq. 
<|A.12fl for a = (a.b) serves to reproduce eq. (|3.6f) . constitutes a posteriori 
verification of the trial solution IjA.llfl . Substituting eqs. I|A.11|) for a = (a,b) 
into eq. ljA.4|) with m — 1 furnishes the first-order total moment, 

i\,r T~T*i , KbB a +K a B b . 
M 1 fvU*t+ +exp, (A.14) 

&a + J^b 



13 



which grows linearly in time at the rate U*. 

Substitution of eqs. I|A.10|) and (| A. 1 1|> for a = (a, b) into eq. i|A.5|) with 
m = 2 yields 



dM 2 
dt 



2 U*U*t + D 



M 



U a K b B a + U b K a B t 
K a + Kb 



exp, 



(A.15) 



wherein D M is given by eq. (|3.9(1 . Use of the latter equation, together with eqs. 
l[A.7j) . f3.6| . (|A.13jl and (|A.14jl . furnishes the particle dispersivity (|3.8j) . 
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Figure Captions 

Figure 1. Comparison between the numerical solution, P (solid line), of the 
microscale equations and the solution, P (dashed line), of the macrotrans- 
port equation for the relatively short time i = 0.5, and for four different 
choices of the equation parameters: (a) K = 1.5, U = —2.0, <t> a = 0.75; 
(b) K = 2.0, U = 0, 4> a = 0.25; (c) K = 0.75, U = 0.5, 4> a = 0.33; (d) 
K = 0.5, U = -0.5, 4> a = 1. 

Figure 2. Comparison between the numerical solution, P (solid line), of the 
microscale equations and the solution, P (dashed line) , of the macrotrans- 
port equation for the relatively long time t = 50, and for four different 
choices of the equation parameters: (a) K = 1.5, U = —2.0, 4> a = 0.75; 
(b) K = 2.0, U = 0, (t> a = 0.25; (c) K = 0.75, U = 0.5, <\> a = 0.33; (d) 
K = 0.5, U = -0.5, 4> a = 1. 
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Dimensionless Converted Position,* Dimensionless Converted Position,* 

(e) (d) 
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